Space-time dimensionality T> as complex 
variable: calculating loop integrals using 
dimensional recurrence relation and analytical 
properties with respect to V. 
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^ ■ Abstract 

We show that dimensional recurrence relation and analytical properties of the loop 
integrals as functions of complex variable T> (space-time dimensionality) provide a 
regular way to derive analytical representations of loop integrals. The representa- 
tions derived have a form of exponentially converging sums. Several examples of the 
. developed technique are given. 
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1 Introduction 



High-order radiative corrections to different amplitudes and cross section be- 
come nowadays more and more important both in QED and QCD. In QED 
this is due to the high precision of the modern spectroscopy experiments, while 
in QCD the radiative corrections play important role due to the strength of 
the interaction. The radiative corrections are expressed in terms of the loop in- 
tegrals, so the possibility to calculate the latter is also very important. Several 
powerful techniques have been developed for the calculation of the multiloop 
integrals. Among them is the IBP (integration-by-part) reduction procedure 
[T1I2] . which allows one to reduce arbitrary loop integral to some finite set of 
master integrals. Owing to IBP reduction, the problem of calculation of loop 
integrals is reduced to the calculation of master integrals. It is important to 
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understand that IBP reduction also helps in the calculation of the master inte- 
grals. Indeed, differentiating the master integral with respect to some external 
invariant or mass, we obtain linear combination of integrals of the same topol- 
ogy or its subtopologies. Applying IBP reduction procedure, we can express 
this combination via the same set of master integrals. Acting in this way, one 
can obtain the system of ordinary differential equations for all master inte- 
grals [HE]- The general solution of this system depends on several constants 
parametrizing the solution of the corresponding homogeneous system. These 
constants can be fixed from the value of the integrals at some specific kine- 
matic parameters, or from their asymptotics in which one invariant becomes 
big, or small, in comparison with others. This method can be applied to the 
integrals with several external parameters. 

However, the integrals with one scale are also important for the applications. 
For them, the method of differential equations does not work and one has to 
rely on other methods. In addition to the direct calculation using, e.g., Mellin- 
Barnes representation, one can try the method of difference equations. One of 
the variants of the difference equation is described by Laporta in Ref. [5]. In 
this method one considers the generalized master integral which is formed by 
raising one of the massive denominators of the master integral of interest to 
arbitrary power x. Using IBP reduction one can obtain the difference equation 
with respect to x. The part, corresponding to the general solution of the ho- 
mogeneous equation is then fixed from the large-x asymptotics. The resulting 
expression for the master integral has the form of factorial series with power- 
like (harmonic) convergence (the general term of the series falls down as k~ a , 
where k is the summation variable, a > 1). This method has been successfully 
applied for several multiloop tadpole integrals, see, e.g., Refs. [6fT] . Due to 
the power-like convergence of sums, the difficulties of obtaining high-precision 
results rapidly grow with the number of digits. The multi-digit results, on the 
other hand, are important for the pslq algorithm [T2], which allows one to 
restore the analytical result from the high-precision numerical one, provided 
the basis of transcedental numbers is known. 

Another variant of the difference equation is the dimensional recurrence re- 
lation connecting the master integrals in different dimensions, suggested by 
Tarasov in Ref. [8]. We remind, that, in contrast to the differential equation, 
where the solution of the homogeneous equation is parametrized by several 
arbitrary constants, in the case of difference equation the corresponding solu- 
tion is parametrized by several arbitrary periodic functions. In Refs. [8|9] the 
large-P asymptotics has been used for fixing these functions. This asymptotics 
has been derived from the explicit parametric form of the integrals. However, 
for multiloop integrals the calculation of this asymptotics appears to be hardly 
accessible. This is why the method of the dimensional recurrence relation was 
mainly applied to the calculation of one-loop integrals. It is worth noting that 
for the case of several external invariants, one may rely on the combination of 
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differential and difference equations methods, see Ref. [TO] . 

Recently, in Ref. [11] the method of the dimensional recurrence relation was 
successfully applied to the calculation of certain four-loop tadpole. In that 
paper, an arbitrary periodic function parametrizing the solution of the ho- 
mogeneous equation was numerically shown to be equal to zero by using the 
Laporta-like difference equation with respect to the power of massive denom- 
inator. However this vanishing of the homogeneous part was rather an excep- 
tion than a rule, and in many known cases the homogeneous part is not zero. 
Needless to say that in these cases the correct choice of the periodic functions 
not only can be hardly justified, but plainly difficult to guess. On the other 
hand, all examples of application of the dimensional recurrence relation show 
that results obtained with the help of this method have a form of exponen- 
tially converging series (the general term of the series falls down as a~ k , where 
k is the summation variable, a > 1). The rapid convergence of these series 
allows one to use the pslq algorithm to express their e-expansion in terms of 
conventional transcedental numbers, like multiple zeta-f unctions. 

In this paper, we describe a method of the calculation of loop integrals based 
on the usage of Tarasov's dimensional recurrence relation. The key point of our 
approach is the use of the analytical properties of the integrals as functions 
of the complex variable T>. These properties allow one to fix the periodic 
functions up to some constants which can be found from the calculation of 
the integral at some definite T>. 

The paper is organized as follows. In the next Section, we introduce notations 
and formulate a general path of finding the master integrals. In Section III 
we briefly analyse the analytical properties of the parametric representation 
needed for our consideration. In Section IV we give several examples of the 
application of the formulated approach. Section V contains discussions of the 
generalization of the formulated approach to the topologies with more than 
one master integral. The last Section gives short conclusion. 



2 Dimensional recurrence relation 



Assume that we are interested in the calculation of the L-loop integral de- 
pending on E linearly independent external momenta p\, . . . ,p^- There are 
N = L(L + l)/2 + LE scalar products depending on the loop momenta U: 

Sik = h-qk, i = l,...,L, k = 1, . . . ,L + E, (1) 
where q x ^ L = h,...,L, Ql+i,...,l+e = Pi,...,e- 
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The loop integral has the form 



j<P) ( n ) = I " " " (2) 



^LV/2 D nx D n 2 p 



n N 
N 



where the scalar functions D a are linear polynomials with respect to s^. The 
functions D a are assumed to be linearly independent and to form a complete 
basis in the sense that any non-zero linear combination of them depends on 
the loop momenta, and any Sik can be expressed in terms of D a . For further 
purposes it is convenient to introduce the operators A a and B a , see Ref. [13] . 
acting as follows 

(A J {7?) ) {n Xl . . . , n N ) = riiJ^ (m, . . . , + 1, . . . , n N ) , 

(BijW) (n 1} ...,n N ) = J (v) (n u . . . , m - 1, . . . , n N ) . (3) 



Relation between the loop integrals in different dimensions has been first in- 
troduced in Ref. [BJ. It was derived using the parametric representation for 
the loop integrals. For the integral given by some graph, the result of Ref. [8] 
may be represented as follows 



J^- 2 Hn) = ^ L j:(A il ...A iL J^) (n), 



(4) 



trees 



where ii,...,«£ numerate the chords of the tree, and /x = ±1 for the eu- 
clidean/pseudoeuclidean case, respectively. It is natural to label this relation 
as a raising one since it expresses any integral in T> — 2 dimensions via the 
integrals in T> dimensions. It may be also convenient to use the lowering di- 
mensional relation which we will obtain now. Our derivation is based on the 
Baikov's approach [T3] in which one passes from the integration over loop mo- 
menta to the integration over scalar products . The master formula can be 
represented as 



dPh . . . d D i L 
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LV/2 



f 



T[(V-E-L + 1)/2,...,(V-E) /2] 

l l+e N [v(i 2 ,...i L , P1 ,..., PE )r- E - L -^ 



j nn. da 



i=l j=i 



[V (Pi,...,Pe)) 



(V-E-l)/2 f ' 

(5) 



where 



V(q 1 ,...,q M ) 



qi ■ qu 



Qi ■ Qm ■ ■ ■ 

is a Gram determinant constructed on the vectors q±, . 



(6) 



qu and / is arbi- 



trary function of scalar products s,,-. For our purposes, of course, we choose / 
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[D^D^ 2 . . . D^]~ . Due to the aforementioned properties of basis {Di, . . . , D^}, 
the function V (Z 2 , . . . Il, Pi, ■ ■ ■ , Pe) has the form of some polynomial of degree 
L + E of Di, . . . , Dn (and of external invariants, which is not important for 
our purposes): 

V(l 2 ,...l L ,p 1 ,...,p E ) = P(D 1 ,D 2 ,...,D N ) (7) 

Replacing T> — > V+2 in Eq. (JHJ), we obtain the lowering dimensional recurrence 
relation 



Comparing Eqs. (j3J) and (jSJ) we see that in order to obtain the former, one has 
to analyse the graph corresponding to the loop integral, while the latter can be 
obtained without any reference to the graph. So, for computer implementation 
the lowering dimensional relation can be more convenient. On the other hand, 
Eq.(jl]) is more "economic" since the right-hand side of it contains integrals 
with indices shifted by at most L, while the right-hand side of Eq. (jSj) may 
contain integrals with indices shifted by L + E. Of course, after IBP reduction, 
both raising and lowering recurrence relations are equivalent. 

In case when the integral in the left-hand side of Eq. (j3j) or Eq. (JSJ) is master 
integral, IBP reduction of the right-hand side gives a difference equation for 
this master integral. If the topology contains only one master integral, the 
right-hand side, after IBP reduction, contains this master integral and master 
integrals of the subtopologies, which we assume to be found by the same 
method, or some other. In case of several master integrals in the topology, 
we can obtain in the same way a system of difference equations for them. 
We will concentrate below on the case of one master integral in the topology 
of interest. The discussion of the case of several- master topology is given in 
Section V. So, after the IBP reduction we obtain from, e.g., raising relation 
the following dimensional recurrence relation 

j(P-2) =C (V)J<P) + R(V), (9) 

where R (V) is non-homogeneous part constructed of the master integrals of 
subtopologies in T> dimensions. The coefficient C (P) is a rational function of 
D which can always be represented as 

C( P) = e n -jr^ 2 ! (10) 

where c is some constant, and a, (hi) are the zeros of the numerator (denomi- 
nator). We can construct the solution of the homogeneous part of the equation 
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as follows: 



T {V) 



cj (V) /£ (V) 



where u (T>) = uj (T> + 2) is arbitrary periodic function of T> and (T>) is 
any specific non-zero solution of the homogeneous equation. One of possible 
choices of £ (V) is 



S- 1 (D) 



_ P/2 nr=ir(q,-xy2) 

n"Li r - P/2) ' 



(12) 



By analogy with the notion of integrating factor in the theory of differential 
equations, we will call S (T>) the summing factor. The choice of the summing 
factor is by no means unique. Given one summing factor, we may multiply it 
by any periodic function of T> to obtain another. In particular, another natural 
choice is 



s- 1 (V) 



' 1 \ n+m 

-1 C 



-p/2 iT7 = irCP/2 + i-y 
nr=i r (v/2 + 1 - a,) 



(13) 



Let us now substitute J {v) = S" 1 (V) g (V) in Eq. (JSJ. After multiplying both 
parts by S (P — 2), we obtain 



g(V-2)=g(V)+r(V) 



(14) 



where r (T>) = R (V) E (T> — 2). Suppose now that we can decompose r (X>) 
into two, r + (T>) and r_ (X>), decreasing fast enough at ±oo, respectively: 



k— >oo 



r(X>) =r+(D-2)+r_(P), r± (V ± 2fc) < a fe , |a| < 1 



(15) 



Note the shift of the argument in the definition of r + (X>). This decomposition 
also determines the corresponding decomposition 



R (V) =R + {V-2) + R_ (V) 



(16) 



such that 



lim 



C(V + 2) R + (V + 2) 



lim 



R + (V) 
C- 1 (V - 2) i?_ (V - 2) 



R- (2?) 



< 1, 

< 1. 



(17) 
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The general solution reads 



r + (V + 2k) 



k=0 k=0 
V 



s + (V, k) = ' + ^-p = cj " ^— -R + (V + 2k) 



n f + i-o* 



m , s 

r (P - 2fc) , fci Si ^ " ^ 



fc+i 



s _ (X>, jfc) = - v 7 = c- 1 -"^ R_ (D - 2k) . 



fe+i 



Note that the specific nonhomogeneous solution does not depend on the par- 
ticular choice of the summing factor. Eq. ffT8l) clearly demonstrates that, in 
order to find J^ v \ we have to fix the periodic function u (T>). This function 
can be conveniently understood as a function of the complex variable 

z = e™ v (19) 

and below we switch the notations u (T>) — > u(z). Understanding u as a 
function of z, we automatically account for the periodicity due to identity 
e «7r(©+2) _ e «7rD_ Q ur [^ ea Q f fixing u) (z) is simple and yet quite efficient as we 
will demonstrate in Section IV. The functions r± in the right-hand side of Eq. 
(1181) are assumed to be known and, in particular, their analytical properties are 
known. The integral in the left-hand side is not known (this is the purpose 
of our investigation), but we can discover some of its analytical properties 
from, e.g., its parametric representation. Note that we do not have to analyse 
the integral in the whole complex plane of T>, but only in an arbitrary stripe 
of width 2, i.e., in S — {T>\ d < KeV ^ d + 2}, which we will call the basic 
stripe. The factor E (T>) is also a known analytical function of T>. This allows 
us to fix to (z) up to some constant which can be found from the value of 
the integral at one given value of T>. There is a number of technical 
conveniencies in this approach. First, as we have said above, the summing 
factor S (T>) can be chosen in different forms and we are free to choose the 
one with the most convenient analytical properties. Second, we can pass from 
using one integral of certain topology as a master to using another integral of 
the same topology as a master. Finally, we are free to choose the convenient 
basic stipe in which we perform the analysis. 

So, summarizing, the suggested path of finding a master integral is the follow- 
ing (we assume we can do IBP reduction): 

(1) Make sure all master integrals in subtopologies are known. If it is not so, 
start from calculating them. 

(2) Pass to a suitable master integral. It is convenient to choose a master 
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integral which is finite in some interval T> G (d, d + 2). For this purpose, 
e.g., increase powers of some massive propagators. 

(3) Construct dimensional recurrence relation for this master integral. Due 
to step 1, the nonhomogeneous part of the recurrence relation is known. 

(4) Find a general solution of this recurrence relation depending on function 
u{z). 

(5) Fix the singularities of this function by analysing the analytical properties 
of the master integrals and summing factor. The convenient choices of 
summing factor and of basic stipe are at your disposal. 

(6) If needed, fix the remaining constants from the value of the integral at 
some space-time dimension T>. 

The step 5 is the key point of our approach. We will present the applica- 
tions of the formulated approach on several examples, but first we would like 
to demonstrate how some of the analytical properties of the integral can be 
determined from its parametric representation. 



3 Parametric representation 



Let us consider the parametric representation of some L-loop integral in Eu- 
clidean space with / internal lines (see, e.g., Ref. [15J): 



r \0 (x)]' DL/2 ~ I 

= r (J - LP/2) J dx x . . . d Xl 5 (1 - Zxi) rpT \-iV(l+i)/2-i ■ (2°) 



[P(x))' 



The polynomials Q (x) and P (x) are determined by the graph. For our con- 
sideration it is important only that both these functions are nonnegative in 
the whole integration region, 



Q(x)^0, P(x)^0. 



Suppose now that the parametric representation converges for all T> in some 
interval (d\, <i 2 )-Then it is easy to see that Eq. ([20} determines J (c) holo- 
morphic function on the whole stripe S = {V\ ReD £ (di,^)}- Moreover, 
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as ImD tends to ±00 while ReL> = g?G (d\, <Lj) is fixed, we can estimate 

[q i^r 12 - 1 



= \T{I-LV/2)\ j dx 1 ...dx I 8{l-Y.x i ) 
< |r (I - LP/2) I J dx 1 ... dxrf (1 - J>i) 
= |r (/ - LP/2) I J dxx... d Xl 5 (1 - Ex,) 



[P (x)]^ 1 )/ 2 - 7 

[g(^,p,m)] M/2 - 7 



[P (a;)] I,(L+1)/2 - / 

[Q^P .m)]^- 7 
[P (*)] 

|r(/-LP/2)| j(d) congt x e _-^ |Imp| ,-l/ 2 -LRe C /2 



d(L+l)/2-I 



r (J - Ld/2) 



f21^ 



These properties of the parametric representation will be used in the examples 
in the next Section. 



4 Examples 



In this Section we give several examples of the application of our method to 
different loop integrals. Except for the last example, our consideration is re- 
stricted by the one-scale integrals, since several-scale integrals can be treated 
by the differential equation method. In all cases, the IBP reduction has been 
performed using the Mathematica program, based on Ref. [13]. All integrals 
considered in this paper are known due to some other techniques. Our rep- 
resentations for one-scale cases, however, differ from the known ones, being 
exponentially convergent sums. 



4-1 Three-loop sunrise tadpole: basics of the technique 



Let us consider the three-loop integral 



{V) 




d v kd v ld v r 
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[22) 



This integral has been investigated in Refs. [T6]fT7] and the result is expressed 
in terms of 3P2 [. . . 1 1] . We will follow the path formulated in the end of Section 
II. 

(1) There is one master integral in subtopologies: 



J la 




r 3 (i -v/2) 



(23) 



9 



(2) The integral j[ is a holomorphic function in the stripe S = {T>\ ReD £ [—2, 0)} 
which can be easily understood from the parametric representation: 



4 V) = r (4 - 3D/2) / dxtdztdxsdxJjl-Zxi) ^ 

It easy to realise that any Euclidean integral with all internal propagators 
massive is holomorphic in the whole half-plane Re£> < 0. 

(3) The dimensional recurrence reads 

{v _ 2) _ (3P - 10)3 (V - 2) (v) _ (llP-38)(P-2) 3 {v) 

1 ~ 128(25-4) 1 64(£>-4) la 1 J 

where (a) n = a (a + 1) . . . (a + n — 1) denotes the Pochhammer symbol. 

(4) We choose the summing factor as 

= r(3/2-P/2)r(3-3P/2) 

v ; r(2-p/2) v ; 

Note that the chosen summing factor has neither poles nor zeros in S. 
Representing = (£>) g (V), we obtain 

g(V-2) = g(V) + r(V) 

r (V) = (l^-38)r 4 (2-P/2) 
1 J 4 p r(5/2-P/2)r(6-3P/2) 1 ; 

The nonhomogeneous part r (T> — 2k) falls down as (16/27) fc when k goes 
to 00, so we can write the solution in the form 

1 - (llP-16-22fc)r 4 (l + A;-P/2) , 

+ 16r (2 - P/2) (3/2 - 2?/2) fc (3 - 3P/2) 3fe 1 j 

(5) Let us consider analytical properties of functions entering Eq. (l28l) in S. 

(T>) 

The integral J} is a holomorphic function in S. When ImX? — > ±00, it 
falls down at least as 

jf )| < |r (4 - 3P/2)| ~ e"^ |ImP| 7/2 - 3Re ^ 2 (29) 

which can be clearly seen from the parametric representation (j2"4"j) . The 
summing factor is also holomorphic function in S and 

|S(P)|<e^|ImP|- 2+3RcC/2 (30) 
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in the limit ImP-> ±00. Finally, the specific solution t (V) of the non- 
homogeneous equation is also a holomorphic function in S. and 



37r|Im T> 



l*P)l^e —\lmV\ a (31) 

Note that the limits ImD — > ± 00 correspond to the limits z — ► 0, 00. So, 
from Eqs. (1251), (125]) , (1501). and (13"T]) we conclude that 

u{z) z ^°° |ImX>r, (32) 

where 1/ is some real number not essential for our consideration. Taking 
into account that 

llmPT llmPl" 
hm — — tj — = hm — = 33) 

js—00 \ Z \ a z^O \ z \- a V ' 

for any v and any a > 0, we conclude that (#) is holomorphic function in 
the extended complex plane of z, except, may be z = and z = 00 and 
growing slower than any positive(negative) power of \z\ when z tends 
to infinity (zero). This is sufficient to claim that to (z) is holomorphic 
function in the extended complex plane of z, being therefore a constant. 
(6) We can fix this constant by the condition J : f D_0 ' ) = 1, which, e.g., follows 
from the parametric representation. We finally obtain 

T (V) P r(3/2-P/2,3-3P/2) 1 - (l + llfe/8) (k\)\ k 
1 " T(2-V/2) 0Ffr o (3/2) fc (3) 3fc 

1 " (IIP - 16 - 22k) T 4 [k + 1 - V/2] k 

16T(2-P/2)jfe (3/2 - V/2) k (3 - 3P/2) 3fc ' 1 1 

It is easy to understand that general terms of both sums in Eq. (|34l) 
fall down as (j^) thus providing fast convergence of sums which can 
be conveniently determined as l/log 10 (27/16) ~ 4.4 terms per decimal 
digit. We may have expressed our result f[3~4"|) in terms of p F q . . . ||| , 
but we consider the sums to be more treatable than these functions. In 
particular, in order to obtain the e-expansion of 4 4 ~ 2e \ we can expand 
under the sum sign. It is easy to check that, up to a O (e°), this expansion 
is given by the first term of the second sum: 

Jt 2e) = --1 w Ue w 3 ^7 ^ M + O (e) (35) 

1 3(2e-l)(3e-l)(3e-2)(e-l) U W V ; 

In fact, the next 5 terms of e-expansion ofjj 4 - 2 ^ can be found in Refs. [HUTT 
In order to demonstrate the fast convergence of the obtained representa- 
tion we present here the 40-digit values of coefficients of jf ^ expansion 
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up to (e 6 ) obtained in less than 2 minutes on Mathematical 

^4-26) = +2 .00000000000000000000000000000000000000 x e~ 3 
+ 4.203372677257469503027594126172252080413 x e~ 2 
+ 12.15744322207889158433285265524001769362 x e" 1 
+ 10.34393554171240616956324844514078698433 x e° 
+ 22.40040337934740516273278836060402645950 x e 

- 192.9905185401601934396793596956955193083 x e 2 

- 298.2989210535372236865167423656517391923 x e 3 

- 3327.931223253862483955917547688685145911 x e 4 

- 5027.397388051400645438297589076204952450 x e 5 

- 37321.16865675250679268222655965209414217 x e 6 

+ O (e 7 ) (36) 
4-2 Four-loop tadpole: Dealing with masless internal lines. 



Let us now consider the following four-loop tadpole integral: 



Jo 



d D kd D ld D rd D p 



n 2V k 2 l 2 r 2 



(k + pf + 1 (l+pf + 1 {r+pf + l 



(37) 



The explicit expression for this integral and expansions near D = 4 and T> = 3 
can be found in Ref . . In this example we will demonstrate how to deal with 
the massless propagators. 

(1) There is one master integral in subtopologies: 



J. 



(V) 



r (4 - 3D/2) r 2 (3 - v) r 2 (i - v/2) r (v/2 - 1) 

2a - - r (6 - 2V) 

(38) 

(2) Due to the presence of massless propagators, there is an infrared diver- 

(T>) 

gence at V = 2 in . It is also easy to see that there is an ultraviolet 
divergence at T> = 3. Thus, the integral is a holomorphic function only 
in the stripe S = {T>\ ReP 6 (2, 3)}. This is not sufficient for our pur- 
poses and we instead consider another integral, with massive propagators 
squared: 



J. 



(V) 



d v kd v ld v rd v p 



7T 



(k + p) 2 + l\ [(l+p) 2 + l\ [(r + pf + l 

(39 



12 



Due to IBP identities, the original integral J 2 can be expressed via J 2 
and J^a as follows: 



(p) _ 3(3X?-11)(3X?-10) ~ {V ) 
2 4(2?-4)(2?-3) 3 (22?-7) 2 



3(1? - 2) (32? - 8) (132? 2 - 882? + 148 
128(2? -3) 2 (22?-7) 2 



^ (40) 



~f:D"i 

Thus, we can restrict ourselves to the calculation of the integral .This 
integral is obviously a holomorphic function of 2? when Re 2? e (2,4|). 
This is sufficient for our purposes. So we choose the basic stripe as 
S = {2?| Re 2? G (2,4]}. Let us make one observation which essentially 
simplifies the following consideration. The two terms in the right-hand 
side of Eq. fiHj) have poles of the third order at T> = 3, while the integral 
J 2 clearly has only first-order pole at 2? = 3. Thus, from Eq. (j40j) . we 
can extract two terms of expansion of near T> = 3: 

Jt 2e) = y + (11 - 4 7 - 8 In 2) + O (e) (41) 

(3) The dimensional recurrence for reads 

P>) 4(22?-8) 4 (2?-3) 3 (2?-l) 2 ~ p+2) 
J2 = 3(31? -11)5 h 

r (2?) = - nU2V iV 7) 5 (3299I?6 " 59493p5 + 4440981,4 

- 17561642? 3 + 38798002? 2 - 45402241? + 2198784) J^ +2) (42) 

(4) It is easy to establish from the explicit form of known integral J 2a that 
in this case R- (I?) = 0. The summing factor obeys the equation 

£ (2? + 2) 4(21?- 8)4(1? -3) 3 (2?-l) 2 
E(2?) 3(31?- 11)5 1 ' 

As we noticed above, we can use different forms of the summing factor 
and we choose 



sin (~V - sin f f 2? - \ ) sin 2 ( *f 

E (2?) = ^ ^ ^ (44) 

2^sin(f2?-f)sin(f2?-f) 



r i 


(7-f) 


!!-(?-§) 


1 r(2? - 1) 


r(9 - 2i?)r(5 - i?)r (f - ^) 



Our choice of the summing factor needs some explanations. When choos- 
ing the summing factor, we used the following "rules of thumb": 
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(a) " No undetermined singularities in basic stripe" : The function E (T>) 
is chosen in such a way that E (T>) J 2 has only known singularities 
in S. We did not eliminate the simple pole at V = 3 because the 

(■T*\\ 

singular part of E (V) J 2 near V = 3 is known due to Eq.pTj). 

(b) " Good behaviour at infinity" : We require that the function E (V) J 2 V ^ 
obey the following constraint 

E(P)£ ) Img ^ OOQ 

exp [tt \V\] v ; 

when Re T> G (2,4]. It is very easy to establish from the Feynman 
parametrization, that 

4 £ r[9-2X>] llmPl", 

so we can easily obtain the desired estimate for E (T>) J 2 . We also 
require r + (V) = E (V) R (V) to obey the same constraint 



exp [tt \V\] 
for Re V > 2. 

(c) " Pole minimization" : The function r + (T>) = E (X>) i? (T>) has some 
poles. We try to minimize their orders and/or number in the region 
Re T> > 2 by adding some factors of the form sin ((P — d) tt/2). This 
requirement explains appearance of sin functions in the numerator 
in Eq. (1441) . At some point we were not able to add sin factors due 
to the condition b. 
The pole structure of 



3tt 5 / 2 2 2V ~ 15 Vcos\ 




1 r(7 - 3x»)r (4 - 


<n 


V l\ 2 
2 2) 


sin(vrD) tan - 




1 tan (f + |) 


1 r(8 -2£>) 2 r 


(f + l) 



r + (V) 

x (3299X? 6 - 59493£> 5 + 444098P 4 - 1756164D 3 

+ 3879800P 2 - 4540224P + 2198784) (47) 

is the following. There are 5 series of simple poles in 

V = V l (k) , V 2 (k) , V 3 (k) , V 4 (k) , T> 5 (k) , k = 0, 1, . . . 
T>\ (k) = 2§ + 2A;, V 2 (k) = 2§ + 2A;, 
£> 3 (fc) = 3I + 2ib, £> 4 (Jfe) = 3§ + 2Jfe, 

£> 5 (fc) = 4 + 2A:, (48) 

and also a single pole at V = 3. This pole structure is depicted in FigJTJ 
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60 
40 
20 



n 



n 




Fig. 1. The function r + (D). 
The general solution of the recurrence relation reads 



J<p = E- 1 (V) 



oo(z) + J2 r+ (V + 2k) 

k=0 



uj e 



AttT> 



(49) 



is the arbitrary periodic function of T> we need to 



where uj (z 
fix. 

(5) Let us now analyse the singularities of uj (z). Rewriting Eq. 



as 



uj (z) = E (V) J 2 (C) - ^ r+ (X> + 2fc) 

fc=0 



(50) 



we see that singularities of this function are determined by the singulari- 
ties of the terms in the right-hand side. Due to our choice of the summing 
factor (see Hall , the first term has only known singularities in S. Namely, 
it has simple pole at V = 3 with the residue, determined from Eq. (14T1) : 



£ (3 — 2e) 4 



(3-26) 



3tt 5 / 2 
256e 



+ O (e°) 



(51) 



From the pole structure of r + (X>) we conclude that the first term of the 
sum in the right-hand side of Eq. fl50|) also has a pole at T> = 3, while the 
other terms of the sum have no singularities in this point. Moreover, 



3tt 5 / 2 / \ 



(52) 



so the singularities of the right-hand side of Eq. (15171) at V = 3 cancel. 
Thus, uj (z) has no singularities in z = exp [iit x 3] = — 1. The five series 
of poles (JHD give rise to the poles of uj (z) in 



e iiz/3 e «27r/3 e -i2n/3 g-iTr/3 ^ 



(53) 
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From the constraints (145]) and (1461) we conclude that uj (z) is regular at 
z = 0, oo. Thus, the general form of uj (z), mimicking the structure of the 
poles in the right-hand side of Eq. (1501) , is the following 



uj (z) = a + 



a i 



«2 



piir/3 



pi2ir/3 



+ 



«3 



«5 



g— i2ir/3 2; g — g J 



(54) 



It is convenient to rewrite this form in terms of cotangents. Using the 
identity 



cot ( - (P 



Pn 



. z + e 



(55) 



we represent 



uj (z) = b + 6i cot § (P- 2|) + 6 2 cot § (P - 2 



+ 6 3 cot I (p - 3|) + 64 cot \{V- 3f ) + 65 cot I (£> - 4) (56) 

The constants &!_ 5 are fixed by the singularities of the sum in the right- 
hand side of Eq. (l50|) : 



b; 



71 



Y Res r+(V) 
2 to v=^) 



(57) 



where T>i (k) are determined in Eq. fl48p . They can be expressed as linear 

Using a guess bi = fiH 7 ^ 2 with fi being 



combinations of p F q 



39 



rationals, we were able to establish the simple form of these coefficients: 

3tt 7 / 2 



bi = -b 2 = -63 = 64 = 



4^ 



32 



(5f 



where = denotes the equality checked numerically with at least 10 3 digits. 
(6) The last coefficient b in uj (z) can be shown to be equal to zero by cosider- 
ing the limit of the right-hand side of Eq. (j50p at T> — > 3 and using Eq. 
f l4T|) . Note that in this limit only the first term in the sum contribute. 
Thus, we obtain finally 



J. 



(V) 



XT 1 (X>) 



3tt 7 / 2 
32 



cot § ( P 



Z 3 



cot f ( £> 



2^ 

Z 3 



cot f (P - 3^) + cot § (P - 3|) - 4 cot § (P - 4) 



+ £ r+ (P + 2A;) 



fc=0 



(59) 



where £ 1 (P) and r + (P) are determined in Eqs. f j4*4l) and ( 1471) . respec- 



3 9 



. Thus, the rate 



tively. Again, we can express the sum via p F q . . . 
of convergence in representation (|59l) is l/log 10 (3 9 /2 14 ) w 12.5 terms 
per decimal digit. We have checked numerically the agreement of our 
representation with the known expansions around P = 4 and P = 3. 
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4-3 Two-loop massless onshell nonplanar vertex: Two- directed summation 



This example demonstrates how the method works in the case when both R + 
and i?_ are not zero. Let us consider the integral 

jf) % f d y Vr (60 ) 

«✓'-' J n v l 2 r 2 (I — p) (r — q) (I — r — p) (r — I — q) 

with 

p 2 = q 2 = 0, ( p + q ) 2 = -l. (61) 

In contrast to the previous examples, here all scalar products are pseudoeu- 
clidean, e.g., I 2 = Iq — l 2 . This integral in arbitrary dimension has been cal- 
culated in Ref. [H] and the result is expressed via q +\F q [. . . |1] with q — 3, 4. 
Let us again apply our technique. 

(1) There are two master integrals in the subtopologies: 

j£ ] = = - ^ v ^ J (62) 
4V = = — : ^ \ / ,.).. n \ (m> 

(2) The integral has collinear divergences at T> = 4 and ultraviolet divergence 
at D — 6. So, we choose the basic stripe as S 1 = {D| ReD G (4,6]}. 
Thorough analysis of Feynman parametrization justifies that J3 is finite 
when ReD G (4,6). We also note that at T> = 6 the integral has simple 
pole. 

(3) The dimensional recurrence reads 

4 V) = -f^^f +2) + R + (V) + 2) 

2 (43D 4 - 478D 3 + 1963D 2 - 3530D + 2352) < v) 

R + \ V ) = 7^ TTT^ 773 J 3a 



r(3-p)r(f - i 1 


3 

I 




2 r(D-3) 2 


r(f-s) 
T(4 - V)T (2 - f ) 


ri 


r(z>-2)ri 


(f-4 


) 



(D - 3) (V - A)' 
W 2 + 8581: 
(3D - 8) (D - 4) 5 



fl_ (P + 2) = - 2 ' 37P3 " 313P2 + ^ - 752 ' Jg». (64) 



Here we used the lowering recurrence relation, thus, appearance of Jg a £ in 
the right-hand side. It is easy to check, that R± (V) satisfy the conditions 

(4) We choose the summing factor as 

E (D) = 4^(4 - D) sin (|(D - 5)) sin 2 (|(D - 4)) r (v - ^) . (65) 
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Note the quadratic zero at V = 6, which provides that 



j{T>) _ 



hmS(P)jr=0, (66) 



and, in particular, that the function £ (V) J\ is regular at V — 6. The 
general solution reads 



jf ) = s- 1 (Z>) 



w (^) + £ r + ( P + 2A ~ E r - ( V ~ 2k ) 

k=0 k=0 



(67) 



V^sin(7rP)r(f-f) 


r 


If" 


2) 


2 r 




(v - 3)r 


< 3V 
K 2 


-3) 





X 



(43£> 4 - 478£> 3 + 1963£> 2 - 3530£> + 2352) 



7r 2 2 2£>-6 r / 2 _ 3) r [V - — 

r_ (V) = ^ U — ^ ^ (37£> 3 - 535V 2 + 2554V - 4016 

1 ; (5-P)r(f-6) V 

(5) Due to our choice of the summing factor, the function r + (V) does not 
have any singularities in the region KeV > 4. The function r_ (V) has 
simple poles at 

V = V 1 (k) , V 2 (k) , 5, 6 
V\ (k) = 5\- 2k, V 2 (k) = A\ - 2k, A; = 0,1,... (68) 

Taking into account that E (V) is holomorphic in S, we obtain 



to 



(V) = b + h cot | (V - 5|) + b 2 cot f (D - 4| 
+6 3 cot |(X> - 5) + 64 cot ! (£> — 6) , 



2^ 1 " 4 ^ 2 

00 00 

61 = - Y Res r_ (V-2k), b 2 = - V Res r_ (D - 2fc) , 
63 = - Resr_ (P) = 256tt 7/2 , 64 = - Resr_ (V) = 1280vr 7/2 (69) 

Again, making a guess &i >2 = fi :2 7c 7 ^ 2 , supported by the values of 6 3i4 , we 
check with the accuracy as high as 10 3 digits that 

bl I b 2 = -2567T 7 / 2 (70) 

(6) The constant b can be fixed using the condition f)66p and the properties 
of r±: 

r + (6 + 2Jfe) = 0, r_ (6 - 2e) - ^ = O (e) 

We obtain 

00 

6 = £r_(4-2ife) = 0. (71) 

A:=0 
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Therefore, Eqs. (|57|), (p5]h (|70j), and (j7T]) give us the result for j; 
Note that the sums of r± in Eq. ( |67|) can be expressed via 5F4 

and 4F3 



(2?) 





1 16 




1 27. 



I 27 
I 64 



respectively. The arguments of the hypergeometric func- 
tions determine the rate of convergency of the sums: the first converges 
with the rate 4.5 terms per decimal digit, while the second converges 
with the rate 2.7 terms per decimal digit. We have checked numerically 



the agreement of our result for the e-expansion of J3 
sponding result in Ref. [18J. 



(4-2e) 



with the corre- 



4-4 Two-loop two-masses tadpole: Two-scale case 



Let us now examine how the method works for the integrals with several scales. 
Consider the integral 



J 



(23) 




d v ld v r 



it v [I 2 + 1] [r 2 + 1] (I- rf + m 2 



(72) 



Of course, we can obtain the differential equation with respect to m and put 
the initial condition in the point m = 1. However, we can also proceed as 
before. 



^1) Master integrals in subtopologies are 

OO =m v ~*oO = r2 [1-^/2] 



J 4a 



(73) 



(2) The basic stripe can be chosen, e.g., as S = {V\ ReV e [0, 2)}. 

(3) The dimensional recurrence reads 



•J A. 



m 2 (4 — m 2 ) 
(V -2)(V- 3) 



J 



(T>-2) 



+ 



(m 2 - 2 - 2m v - 2 
(V -2)(V- 3) 



r 2 (2-£>/2) (74) 



(4) Let us assume that m < 2 for the moment. The result for m > 2 can be 
obtained by the analytical continuation. We choose the summing factor 

as 

m -v (4 _ m 2 Y Vl2 



E(P) 



After the replacement J\ 



(23) 



r(2 - v) 

g (V) /S (V) we obtain 



(75) 



g(V) = g(V-2)+r(V) 



r{V) 



m 



2 - 2m D - 2 ] T 2 (2 - V/2) 



mP (4-m 2 ) v/2 T(4-V) 



(76) 
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The limit 



v r{V + 2) 
lim ; — - — 

v^-oo r (V) 



1 — j max ( I . ;// ) % I 



(77) 



shows that r (T>) = r_ (T>) and allows us to represent the general solution 
in the form 

4 V) = X- 1 (V)u(z)+t(V) 



t(V) 



r (i - v/2) r (2 - v/2) 



2 (3 - V) 



k=0 \ ^ / 



(5/2 - Vj2\ 



r (i - x>/2) r (2 - x>/2) 



2 (3 - V) 



(m 2 - 2) 2 Fi 



1,2-23/2 
5/2-23/2 



m 



m 
T 



-2m D ~ 2 2 F 1 



1,2-75/2 
5/2-23/2 



(5) The functions S (£>) jj 13 '* and S (P) t (P) are holomorphic in S 1 , They also 
behave well when ImD — > ±oo, so uj (z) = const. 

(6) Fixing the constant by the condition J4 = m~ 2 , we obtain 



m 
T 



uj (z) 



m 



f(0) 



4tt# (2 - m 2 ) 
m 3 (4 — m 2 ) 3 ^ 2 



(75 



One can check that unusual 9 (2 — m 2 ) in this formula cancels the discon- 



tinuity of (m 2 — 2) 2 F\ 



1,2-23/2 
5/2-23/2 



m 2_ ^ 



in i (Z>) . Using the properties 



of the hypergeometric function, we rewrite the final result as 

r(i-f)r(2-f) 



j 





1,3-23 


2 " 

m 




5/2-23/2 


4 _ 



X>-3 

which coincides with the one obtained in Ref. H~9l. 



rrP- 2 TP 

m 2-^1 



1,2-23/2 


2 " 

m 


) 


3/2 


4 





5 Topologies with several master integrals 



We would like to discuss now the applicability of the described method for 
the topologies containing more then one master integral. In this case the di- 
mensional recurrence relation (Q should be understood as the vector equation 
with j( v \ R (V) being the column vectors and C (T>) the matrix. As it follows 
from the theory of difference equations (see, e.g., Ref. [5]), the solution of the 
homogeneous part of ([9]) can be sought for in the form of factorial series and 
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the general form is 

4 V) — F (V) u (z) , (79) 

where F (T>) is the matrix of fundamental solutions and u (z) = to (e mV ^ is 
a column vector of arbitrary periodic functions of T> (with period 2). So, to 
proceed in the same way as before, we have to define the summing factor as 
£ (X>) = F~ x (T>). While the elements of the matrix F (T>) have the form of 
factorial series and their analytical properties can be analysed in the same 
way as those of the nonhomogeneous terms in our consideration, this is not 
so for the inverse matrix F^ 1 (X?). To overcome this difficulty, we can use 
the following transformation. The matrix F (V) satisfies the following matrix 
equation 

F (V - 2) = C (V) F (V) (80) 
Taking the determinant of both sides, we have the equation 



f(V-2) = c(V)f(V), (81) 

where / (T>) = det F (V), c (V) = det C (V). The coefficient c (V) is a rational 
function of T> and can be represented as 



Ui(ai-V/2) . , 

c (V) = — 82 

1 ; Yljibj-V/2) y } 



Therefore, the general form of / (T>) is (cf. with Eq. (JT2 



f ^^ v,2 m^^r iz) (83) 

Now we may use the formula 

F- l (V)= f^WFp), (84) 

where F {T>) is the adjugate of F (T>). This allows us to choose the summing 
factor in the form 



1 ; nr=ir(a 4 -p/2) K J K } 

This formula does not contain the sums in the denominator which allows one 
to perform the analysis of its analytical properties in the same manner as we 
did before. Again, the choice of the summing factor is not unique due to the 
possibility of its multiplication from the left by any periodic matrix Q (z) . The 
examples of application of this technique to the topologies with several master 
integrals will be given elsewhere. 
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6 Conclusions 



We have presented here the approach to the calculation of the master inte- 
grals based on Tarasov's dimensional recurrence relation and on the analytical 
properties of the loop integrals as functions of the complex variable T>. The 
results obtained within this approach have the form of exponentially converg- 
ing series, in contrast to those obtained with many other methods having only 
power-like convergence. Fast convergence of the sums allow one to apply the 
pslq algorithm in order to obtain the representation of the e-expansion of the 
integral in terms of the conventional transcendental numbers, like multiple 
^-functions. We have limited the presentation to the cases in which master 
integrals of subtopologies are expressed in terms of T-functions. In a more 
general case, when the master integrals of subtopologies are expressed via hy- 
pergeometric sums and/or the main topology contains more than one master 
integral, the analysis of the analytic properties, though being more compli- 
cated, can be performed in the same way. Presumably, this approach can be 
applied to the calculation of various master integrals with one scale, including, 
in particular, the three-loop onshell vertex integrals not yet calculated. 
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